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Abstract 

A rigorous  analysis  of  diffraction  by  rectangularly  pro- 
filed, dielectric  gratings  is  presented.  An  analytical  solu- 
tion is  obtained  for  the  total  electromagnetic  field.  The 
approach  presented  can  readily  be  extended  to  computa- 
tions for  metallic  gratings  and  arbitrarily  shaped  surface- 
relief  gratings. 


1 Introduction 

Because  of  their  many  potential  applications,  two- 
dimensional  (2-D)  dielectric  gratings  are  of  great  interest.  In 
these  applications  the  gratings  typically  assume  functions  of 
beam  combiners,  splitters,  or  shapers,  etc.  [1],  [2]. 

The  design  of  such  gratings  inherently  is  an  inverse 
problem,  in  our  case  involving  Maxwell’s  equations  with 
boundary  conditions.  Genetic  Algorithms  (GA)  have  shown 
promise  of  coping  with  this  kind  of  problem.  In  [3],  an  ex- 
cellent overview  is  provided  on  the  application  of  GA  to  en- 
gineering electromagnetics. 

Genetic  Algorithms  are  stochastic  search  and  optimiza- 
tion techniques  with  many  advantages.  They  tend  to  produce 
globally  optimal  results,  and  the  target  functions  are  not  re- 
quired to  be  continuous.  Another  remarkable  property  is  the 
inherent  parallelism  of  GA  so  that  they  are  ideally  suited  for 
implementation  on  a parallel  machine. 

For  a CAD  tool  to  be  useful  in  the  design  of  2-D  grat- 
ings based  on  GA,  an  efficient  method  for  analysis  is  needed. 
With  further  extensions,  the  method  presented  here  has  the 
required  properties.  Modified  slightly,  it  has  been  succesfully 
applied  to  analyze  Multi-Quantum-Well  DFB  laser  diodes  in 
a three-dimensional  manner  as  demonstrated  in  [4]. 

In  [5],  a system  of  second-order  differential  equations 
has  been  introduced  to  describe  the  total  electromagnetic 
field  within  a rectangularly  profiled,  dielectric  grating1.  This 
system  can  be  regarded  as  a wave  equation  with  propagation 
constant.  In  this  paper,  an  analytical  solution  to  this  wave 
equation  is  presented.  For  the  sake  of  convenience,  the  wave 
equation  will  be  briefly  derived. 

l...  also  called  binary  dielectric  grating. 


Figure  1:  A binary,  biperiodic  grating  exposed  to  the  field  of 
an  incident  plane  wave. 


2 Theory 

In  the  scattering  of  a homogeneous  plane  wave  E*(r)  = 
E0e_jk  r,  characterized  by  its  polarization  vector  E0  and  its 
wave  number  vector  k,  by  a biperiodically  structured  object, 
as  illustrated  in  Fig.l,  the  total  electromagnetic  field  obeys 
the  Floquet  theorem,  i.e. 

[E(r  - rmn),  H(r  - rmJ]  = e>'kr"*"  [E(r),  H(r)]  (1) 

for  any  r,  where  rmn  = mLxux  -f  nLyuy.  Hence  the  total 
field  can  be  expressed  by 


E(r) 

H(r) 


M N 

= E E 
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Emn  (#) 
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with  am  = k • ux  + m27r/Lx  and  /3n  — k-uy  + n27r/Ly. 
Lx  and  Ly  denote  the  periodicity  in  ux  and  uy  direction, 
respectively.  The  vector  expansion  coefficients  Emn  (z)  and 
H mn(z)  are  the  unknowns  which  must  be  solved. 

Slicing  a scatterer  in  layers,  the  permittivity  e(:r,  y , z)  of 
each  layer  at  constant  z can  be  expanded  into  a 2-D  Fourier 
series 

M N 

e(x,y,z)=  Y,  (3) 

m=— M n=—N 
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In  a homogeneous  layer,  e(x,  y,  z)  simply  equals  £oo*  Note 
that  in  Eqs.  (2)  and  (3),  M and  N approach  infinity. 

Substitution  of  expansions  (2)  and  (3)  into  Maxwell’s 
equations  yields  after  some  manipulations 

0 = ^U(z)  + u^LCU(z)  (4) 

I(z)  = (-jw0L)“1^U(z)  (5) 

where  the  following  vectors 

ur  = [u*[. . . ,E„(m>n), . . .],  Uy  [. . . . . .]]  (6) 

I = [Ux  [.  . . , H,/(m,n)  !■■*]>  ”uy  [*  • • > ?••*]]  (7) 

with  a bijective  index  transformation  v — (2N  4*  l)(m  + 
M) + n + N 4*  1 are  defined.  Further,  denoting  the  Kronecker 
product  with  the  symbol  ®,  the  identity  matrix  by  1,  and 
letting  k0  = u>0\/A diagonal  matrices 


3 Analytical  Solution 

The  system  of  second-order  differential  equations  (4)  can 
easily  be  solved  by  diagonalizing  [jz]2.  If  P contains  the 
right  eigenvectors  of  [7 z]2,  the  propagation  constants  of  the 
waves  in  the  decoupled  domain  are  given  by  the  following 
diagonal  matrix 

[%]  = ±\/FWP  (15) 

A characteristic  admittance  and  impedance  can  be  defined  as 
is  also  done  in  transmission  line  theory: 

Y„  = Z-1  = P-^jWL)-^  [7,]  (16) 

= P-‘(iu0C)P  [7,]-1  (17) 

The  solutions  to  the  wave  equation  (4)  are  exponential  func- 

tions. After  rearranging  terms,  the  analytical  solutions  can 
assume  the  form 


2N+1 

[a]  = i[(k  • u.)l  + y— diag([— JW, ....  M]  ® [WO])]  (8) 

K0  ■L'x 

= j|4(k  • uv)l  + ^diag([l1_1l]®[-Ar)  ...,N])}  (9) 

2M  + 1 

and  the  block  matrices  become 

1 - [a]N~2[a]  -[a]J\f~2[j3]_ 
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where  the  block  Toeplitz-matrix 
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with 
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Eqs.  (4)  and 

(5). 

(13) 


Evidently,  the  diffraction  of  electromagnetic  waves  by 
objects  with  periodic  structure  can  be  modelled  by  coupled 
transmission  lines  [6],  because  they  are  also  described  by  a 
wave  equation  of  exactly  the  same  form  as  given  in  Eqs.  (4) 
and  (5).  It  has  to  be  pointed  out  that  the  definition  of  the  ma- 
trix for  the  refractive  index  squared,  A/-2,  allows  this  compact 
formulation. 

As  in  transmission  line  theory,  the  square  of  a propaga- 
tion constant  in  z-direction  can  be  defined  in  Eq.  (4): 

[7*] 2 = -ufrC  (14) 

Accordingly,  a similar  wave  equation  with  propagation  con- 
stant -cj2CL  can  also  be  deduced  for  the  magnetic  field. 

Note,  that  the  presented  formulation  is  exact , because 
expansions  (2)  and  (3)  for  M,  N oo  are  both  rigorous. 
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= TA(z)T~1 
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where 


and 


T = 
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0 PY0 
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cosh([7z]z)  sinh([7z]z) 
sinh([7z]2)  cosh([7z]2) 


(19) 

(20) 


Matrix  A denotes  the  chain  matrix  for  one  layer.  For  a ho- 
mogeneous layer,  the  eigenvalues  are  explicitly  given  so  that 
the  chain  matrix  simplifies  accordingly. 

The  electromagnetic  field  in  a structure  consisting  of 
cascaded  layers  is  easily  obtained  once  the  chain  matrix  of 
each  layer  has  been  computed.  If  A*  corresponds  to  the  chain 
matrix  of  layer  i in  Fig.  1,  the  electromagnetic  field  in  a plane 
z,  Z3  < z < Z4  is  then  given  by  the  product  of  the  respective 
chain  matrices: 


r u go 

loo 


A 4(z  - Z3)A3(d3)A2(d2)Ai  (di) 


U(z0) 

I(*>) 


(21) 

In  the  same  manner,  the  total  electromagnetic  field  in  arbi- 
trarily biperiodic  structures  can  be  computed  to  any  desired 
accuracy.  This  is  achieved  through  slicing  the  biperiodic  re- 
gion into  layers  which  are  constant  in  z,  so  that  the  above 
solution  can  be  utilized;  the  thinner  the  layers  the  higher  the 
accuracy.  The  overall  chain  matrix  is  again  calculated  as  the 
product  of  the  chain  matrices  of  the  respective  layers. 


4 Numerical  Example 

In  Figs.  2 and  3,  the  magnitude  of  the  ar-component  of  the 
total  electric  field  is  illustrated  for  the  case,  that  an  inci- 
dent Gaussian  beam  at  65GHz  (beam  waists  wox  = woy  = 
10mm)  is  splitted  by  a binary  grating  as  shown  in  Fig.  1.  The 
values  of  the  design  parameters  are  as  follows:  £\  = £4  = 1, 
e2  = S3  = 2.25,  d\  = 8.7mm,  d2  = 5.8mm,  d3  = 15.0mm, 
c/4  = 30.5mm,  Lx  — 7.0mm,  Ly  = 8.3mm,  bx  = 5.6mm, 
by  = 3.3mm.  The  incident  beam  has  only  an  rr-component, 
so  that  no  ^-component  of  the  total  electric  field  exists.  The 
images  have  been  produced  by  using  a more  elaborated  ver- 
sion of  the  analysis  presented  here. 
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Figure  2:  |^| -component  of  the  total  electric  in  planes 

mmm  Und  g 

5 Conclusions 

Since  all  techniques  applied  in  transmission  line  theory  can 
also  be  applied  to  the  diffraction  of  electromagnetic  flei(js  at 
biperiodic  structures,  the  analysis  presented  still  incorporates 
many  imaginable  extensions  of  the  theory.  For  instance,  re- 
flection anc*  transmission  coeffjCjents  in  matrix  form  can  be 
defined  *n  a similar  fashion  as  has  been  shown  for  the  propa- 
gation constant  and  the  characteristic  admittance/impedance. 
By  introducing  a complex  permittivity,  the  analysis  presented 
here  can  also  be  extended  to  metallic  gratings.  Furthermore, 
the  diffraction  of  electromagnetic  beams  can  be  treated,  if 
they  are  periodically  shifted  [7]. 

To  summarize,  an  efficjent  and  compact  analysis 
method  has  been  presented  which  allows  the  analytical  so- 
lutions known  from  transmission  line  theory  to  be  applied 
to  electromagnetic  diffracted  by  objects  with  periodic 
structure  in  one,  two,  and  three  dimensions.  As  indicated, 
further  generalizations  of  the  theory  are  possible. 
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Figure  3:  j ^ | -component  of  the  total  electric  in  planes 

2 = 0 anc*  y — 0 
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